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We propose a novel trilayer 7r-junction that takes advantage of the unconventional s x 2 y 2 = 
cos k x cos k y pairing symmetry which changes sign between electron and hole Fermi pockets in the 
iron-pnictides. In addition, we also present theoretical results for Andreev bound states in thin 
superconductor-normal metal (or insulator)-iron-pnictide junctions. The presence of non-trivial 
in-gap states, which uniquely appear in this unconventional pairing state, is a distinct feature in 
comparison to other singlet pairing states. 



A family of iron-based high-temperature superconduc- 
tors has recently been discovered^ These compounds 
triggered enormous experimental and theoretical inter- 
est. In particular, probing the Cooper pair symmetry 
is critical to understanding the pairing mechanism of 
this new type of superconductors. Theoretically, many 
possible gap pairing symmetries have been proposed for 
iron pnictides, due to the material's multi-orbital nature 
and complex Fermi surfaces (FSs) , with two hole pockets 
around T-point and two electron pockets around M-point 
[see Fig. [TJ (a)]. 

Among all the candidates, the proposal of s-wave pair- 
ing symmetry with relative sign change between hole 
and electron pockets has appealing advantages'' Two 
of us have predicted^ based on a local magnetic ex- 
change coupling J1-J2 model;" an unconventional s- 
wave symmetry with a particular s x 2 y 2 = cos k x cos k y 
form in the reciprocal momentum space. The pre- 
dicted order parameter is consistent with the relative 
values of the gap on the hole and electron FSs reported 
by angle-resolved photo-emission spectroscopy (ARPES) 
experiments^ The proposed symmetry is also consis- 
tent with low temperature-dependent penetration depth 



experiments, 



and partially explains nuclear spin- 



lattice relaxation ratei 13 i 14 However, since most exper- 
iments arc only sensitive to the magnitude of the gap of 
superconducting (SC) order, a direct p/iase-scnsitive ex- 
periment is essential to map out the complete picture of 
the pairing symmetry. So far there is no proposed di- 
rect phase-sensitive experiment for iron pnictides similar 
to the dc superconducting quantum interference device 
(SQUID) interferometer for the cuprates. The difficulty 
arises from the non-trivial phase structure of the order 
parameter in k space. One possible phase-sensitive ex- 
periment is Andreev spectroscopy in the normal metal 
to superconductor (NS) junction. Unfortunately, two re- 



cent experiment: 



,15.16 



give seemingly conflicting results 



and detailed theoretical study shows indistinguishable 
features between an usual s-wave and sign-changed s- 
wave symmetries^ 

In this paper, we theoretically consider two types of 
Josephson junctions which have novel properties uniquely 
associated with a sign-changed s-wave SC order as op- 



posed to other (singlet) pairing symmetries. Specifically, 
the first type of junction we consider is a trilayer SC de- 
vice where the iron-based superconductor is sandwiched 
by two other layered, s-wave superconductors [see in- 
set of Fig. [1] (b)]. With certain chosen FSs of the two 
outside layers which couple stronger, respectively, to the 
hole and electron pockets of the iron pnictide, due to 
momentum conservation, it is shown that sign-changed 
s-wave pairing symmetry uniquely gives rise to a n- 
junction behavio n 18 i 19 The second type is a single-band 
superconductor-normal metal (or insulator)-iron pnictide 
(SNS'/SIS') junction (see Fig. [2]). Based on the similar 
physics of Andreev reflection at the interface between a 
normal metal and a superconductor ) 20 ! 21 we demonstrate 
that, by adopting a minimal two-orbital model 2 ' to in- 
clude the multi-orbital effect and complex FSs, the non- 
trivial phase structure of the sign-changed s-wave sym- 
metry shows up in the profile of the quasi-particle (QP) 
local density of states (LDOS) in the normal region of the 
junction (see Fig. [¥]): the sign-changed s-wave symmetry 
state supports in-gap bound state solutions. 

Novel ir-junction. We propose a composite Josephson 
junction in which 7r-junction behavior can occur based on 
the unusual phase structure of the s x 2 y 2 -wave pairing. A 
7r-junction defines the situation when the Josephson cou- 
pling J between two superconductors becomes real and 
negative (with no spontaneous or explicit time reversal 
symmetry breaking). In other words, the ground state 
energy (GSE) is minimized as the phase difference <E> be- 
tween two superconductors is 'V , in contrast to the case 
of a "0"-junction. The occurrence of the 7r-shift behavior 
can be usually due to magnetic ordering, strong correla- 
tion effects near the tunneling interface^ or non-trivial 
phase structure of the SC order parameter such as d x 2_ y 2 
pairing symmetry^ 

Unlike these common designs, our proposed junction 
[see the inset of Fig.[T](b)] is composed of an iron-pnictidc 
(S m ) sandwiched by a top and a bottom quasi-2D s-wave 
superconductors (St and S&). The interface between any 
two superconductors is an insulating thin film playing 
the role of a tunneling barrier. The key requirement for 
the top and bottom superconducting materials is that 
the Cooper-pair tunneling probability is stronger into the 
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FIG. 1: (Color online) The plots of (a) schematic stacked 
Fermi surfaces from tri-layer superconducting junction and 
(b) the ground state energy deviation (from $ = case) per 
site as a function of relative phase $ with 6 m — 0. The colored 
curves in (a) represent the FSs of the iron-pnictides. 



hole (electron) pockets for the top (bottom) or vice versa. 
This could be engineered to be due to the normal state 
FSs of the top and bottom superconductors. One possible 
way to achieve this condition is to select a small FS and a 
large FS around T point for the top layer and the bottom 
layer, respectively [see Fig. Q](a)], provided the in-plane 
(perpendicular to the tunneling direction) momentum is 
conserved ideally after tunneling. 

A simple mean-field model Hamiltonian for this tri- 
layer junction can be of the form, Hj = H t + H m + 
H b + H T , where H„ = £ k $f [(e„, k - m)^3 + A v e w ^cr + + 
A l; e~' i6 '' 7 c7_]?/>, ) for n = t,b. e l; . k has the form of 
— 2t n (cosk x + cosky) + e v , o± = {o~\ ± ia-i)/2, and ip v 

is the usual Nambu spinor, xjj^ = (c^ k j, c^-k.j.)- The 
difference between the top and bottom SC phases is 
gauge invariant for the whole junction and is set to be 
<3> = 9 t — 9b. H m , the Hamiltonian of the iron pnictide, 
is shown in Eq. ([5|) transformed into momentum space 
with band parameters and non-vanishing A S 2 = Ae t6m 
given in the caption of Fig. [3] The tunneling Hamilto- 
nian, Ht, which connects neighboring layers, takes the 
simple form: H T = E P ,k,r / ^(p)^ fl '/ , m (k), where hr, v 
is a 2 x 4 matrix, 

^ -.9, -g v ) ^ k (L > 

and the spinor ip m = (c\ M , Ci _ k)i , c\ M , c 2 ,-k,i) de- 
scribes the iron pnictide. Note that we have assumed 
that the dispersion along z-axis is irrelevant and negligi- 
ble in quasi-2D materials. 

For demonstration purpose, we choose parameters 
t t = h = l,e t = 4.78, e b = 1.88, g t = g b = 0.01, 
A t = 0.5, A fc = 0.4, A = 0.5, and 6 m = 0. The stacked 
FSs in the first Brillouin zone from each layer is shown 
in Fig. [T] (a). Now, it is easy to diagonalize Hj and the 
ground-state energy for this mean-field Hamiltonian is 
simply the sum of all QP eigen-cnergies below E = 0. 
As presented in Fig. [T] (b), the ground state energy per 



site relative to the energy of $ = decreases as a func- 
tion of $ with its minimum located at 'V . The physics 
of this result can be easily captured by the perturbation 
calculations of GSE, which give, up to second order of g^, 
—Jtm cos(#t— 6 m ) — J m b cos(9 m —9b), where the Josephson 
couplings, J tm > 0, J m b < (the sign difference between 
the Josephson couplings is due to ir phase difference be- 
tween electron-like and hole-like FSs), have similar form 
of the textbook derivation^ It is obvious to see that the 
global minimum reaches at \9t — 9b\ = 7r. The overall 
junction shows 'V behavior. 

Some comments on the experimental realization are in 
order. First, "large" or "small" FS is meaningful only 
when the lattice constant a is equal or comparable to the 
iron pnictides, where the nearest-neighbor Fe-Fe distance 
is around 2.85A. Second, to make the tunneling processes 
reasonably dominated by in-plane momentum conserva- 
tion, quasi-2D, s-wave SC materials should be used for 
the top and bottom layers due to their irrelevant disper- 
sion along z direction and the epitaxial growing technique 
may be useful for making the coherent-tunnel interfaces. 
Based on these considerations, some plausible candidates 
for the large FS are, for instance, MgB 2 from its ir band 
with a ~ 3A or thin film of Beryllium with a ~ 2.3A; for 
the small FS, it could be 2H-NbSe 2 , where a ~ 3.45A^4 
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FIG. 2: (Color online) A schematic plot of the SNS' junction. 

SNS' /SIS' junctions. A further feature of the s x 2 y 2- 
SC is revealed by considering a Josephson junction which 
connects, on one side, a single-band s-wave superconduc- 
tor, through a normal metal to, on the other side, an iron- 
based superconductor^ We assume its QP spectrum is 
well approximated by a BCS type mean-field Hamilto- 
nian subject to an inhomogeneous pairing field along the 
tunneling direction. Ignoring the z-axis for simplicity, 
the 2D model Hamiltonian of this junction reads, 

H = G(-x-^)H s + G(^-\x\)H N + e(x-±)H s , 
+ H T , (2) 
Hs = -ts ^2 c l,/ r '»"i ii E c !><' c ^ 

<rr'),cr r,a 

+ 5^(A i 4 iT 4 4 + h.c.) ) (3) 

r 

H N = -t N J2 fUfr',. + (U-fl)Y, fL.fr,*, (4) 
(rr'),cr r,<r 

where 0(x) is the Heaviside function, the shifted chem- 
ical potential, p, = fj, — e, guarantees a partially filled 



band, and U denotes barrier potential. To mimic the 
iron pnictide FS, we adopt a two-orbital exchange cou- 
pling model. 2 This leads to a somewhat complicated form 
of Hs> = Hq + Ha , where we separate it into the band 
structure and pairing field parts, 
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where a = 1,2 correspond to and d yz orbitals, 
respectively, and the hopping parameters are given in 
Fig- [HI ((rr'}} represents a next nearest- neighbor pair 
and <fi rr > = 1(— 1) as r — r' = ±x(y). All pairing fields 
considered in TL are set to be real^S- Ao,A s i,A S 2, and 
A^ correspond to intra-orbital on-site, cosk x + cosk y , 
cos k x cos k y , and cos fca, — cos fcj, pairing strength, respec- 
tively, and inter-orbital pairing is ignored due to its small 
contribution as discussed in Ref. 0- Finally, Ht describes 
the tunneling amplitudes across two interfaces around 
±d/2. It can be written as 



H 



t = .9sE c ri.,<T/> 



r L ,ah' L ,<r+gS 



^2 c i,r R ,<rfr'R,<r + h.C., (6) 



where r^,r'i (r^,r'fi) are understood to be coordinates 
across the left (right) interface. 

Before proceeding to compute the QP spectrum and 
the corresponding LDOS, it is important to realize that TL 



is particle-hole symmetric under c 



r. 



(for all fermion operators), and hence its spectrum should 
be symmetric with respect to zero-energy. Taking ad- 
vantage of the translational symmetry transverse to the 
tunneling direction x, TL can be further decomposed 
into a sum of ID Hamiltonians by partially Fourier- 
transforming along (100) surface (y direction). As a con- 
sequence, the whole system is mapped onto an ID effec- 
tive lattice in the form of TL = ^2 k Hio(k y ). Basically, 
this transformation results in effective chemical potential 
and pairing fields with fcy-dcpcndcncc. 

We diagonalize the model Hamiltonian H 1 ]j(k y ) for 
— 7r < k y < 7T on the Nambu basis 'J, where in numerical 
calculations we take the total number of the ID lattice 
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FIG. 3: (Color online) QP spectrum in the SNS' junction 
involving the iron-pnictide with various pairing symmetries. 
The red line represents Andreev bound states in the case of 
s x 2 y 2-wctve pairing. The used parameters are U = 0, fb = —2, 
A s = 0.5, g s = 93' = 1, ti = -1, t 2 = 1.3, t 3 = -0.85, 
ti = —0.85, /i = 1.58, n m = 1, and the pairing strength is, 
respectively, (a) A = 0.5 (n=201), (b) A s2 = 0.5 (n=201), 
(c)A,i = 0.5 (n=201), and (d) A d = 0.5 (n=601) in energy 
units, t s = 1. 



sites n sufficiently large with open boundary conditions; 
the normal metal is always set in the middle of system 
and it is 7i m -site wide. This length is always much less 
than the SC coherence length. In Fig. [3J we show numer- 
ical results of the zero temperature QP spectrum as a 
function of k y for various pairing symmetries of the iron- 
pnictides at \i = 1.58 (electron doped). In addition, to 
visualize the Andreev bound states, in Fig.[2]we also com- 
pute each corresponding QP-LDOS as a function of po- 
sition and energy, D(x,lu) = J2k y ,i k y )\ 2 S(ui - e ; ), 
where i denotes the zth eigenfunction. The magnitude of 
the bulk SC order parameters on both sides of the junc- 
tion is taken to be the same. The ratio of the gap to the 
half band-width of the spectrum is around 0.02 and the 
coherence length is estimated as £ ~ i>f/A s ~ 4a. 

As clearly seen in Fig. [31 the presence of the in-gap 
Andreev bound states with significant weight in different 
k y channels for the extended s-wave cos fc x cos /c y (s x 2 y 2) 
pairing symmetry is a sharp feature distinguishing it from 
other pairing symmetries. Although in addition to (gray- 
color filled) continuum states there are discrete energy 
levels for s, s x 2 +y 2, and d x 2_ y 2 pairing symmetries, they 
are either near the maximum gap edge (s) or only ap- 
pear in certain range of k y ( for s x 2 +y 2, d x 2_ y 2). Espe- 
cially for the latter case, due to the presence of nodal 
points on the electron or the hole pockets, the contribu- 
tion from scattering-state can easily overwhelm that from 
the bound states and can lead to qualitatively different 
QP-LDOS from the case of s^^-wave, where a sharp 
peak appears at the positive subgap energy. 

Furthermore, two observations deserve mentioning. 
First, the features shown in Figs. [3] and 0] do not change 
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FIG. 4: (Color online) QP-LDOS at x = in the SNS' junc- 
tion involving iron-pnictide with various pairing symmetries. 
The used parameters are the same as in drawing Fig. [3] except 
that n = 121 for (a)-(c) and n = 161 for (d). 



much for different doping levels as long as the doping 
concentration is not large enough so that the Fermi sur- 
faces pass the nodal line of cosk x cosk y , i.e., k x = ±tt/2 
and k y = ±7r/2. Second, if the barrier potential U is 
greater than the difference between /j, and the band bot- 
tom, the normal region becomes insulating. This moves 
the subgap peak in the s x 2 y 2 LDOS closer to the gap edge 
without destroying it. 

Can we understand the presence of such non-trivial 
bound states for the s^a^a-wave pairing in a simple way? 
A physical insight for this junction involving such an 
unconventional symmetry can be obtained by treating 
the bands at the electron and hole pockets in the iron- 
pnictide as independent of each other j^i Consequently, a 
simple description for the junction based on Bogoliubov- 
de Gennes (BdG) equations reads: 

[H ,x<T3 + A x (r)a+ + A^(r)a_]^(r) = e A V(r), (7) 

where A — 1,2 are the band indices, and <7j(« = 1,2,3) 
are the Pauli matrices with a± = (<xi ± icr-z) /2 acting 

in the Nambu space, V = [SA v r ),tu(r)] t . Ho,\ con- 
tains the band information of non-interacting electrons 



and A^(r) corresponds to the pairing field in the same 
band A. Along the tunneling direction x, the inhomogc- 
neous A x (r) is modeled by A\(x) — A s e l6s as x < and 



A x (x) 



i6x 



as x > 0. For simplicity, we set 6 S = 



hereafter and keep in mind that for the sign-changed s- 
wave symmetry, 9\ =82 + n. For electrons near FSs, it is 
valid to linearize BdG equations within WKJB approxi- 
mation, ip ~ e 4kj? r 0(r), <j>(r) = (u(r),u(r))* and then the 
BdG equations are now reduced to the form of ID Dirac 
equation if we further take the advantage of translational 
symmetry in transverse direction, 

[-iv Fx d x <j 3 + Ax(x)a + + A* x (x)(T-\(j)(x) = t\4>{x), (8) 

After straightforward calculations with trial bound state 
solutions, u(y) ~ u (vo)e~ 1±x , as studied in Ref. H3, 
the discrete energy level within the gap is given by 
E = ±(A S A X sin e x )/ x /A 2 s + A 2 -2 cos 8 X A S A X , pro- 
vided cos6 x < min(AA/A s , A s /A x ). The pair of solu- 
tions with eigenvalues symmetric with respect to zero- 
energy follows from the particle-hole symmetry of the 
BdG equations. It is clear to see that when superconduc- 
tors on both sides of the junction are in phase (6\ = 0) 
no bound state solution is found, while when they are 
out of phase (9\ = it) there are doubly-degenerate zero 
modes i2I 

The significance of this simple result is that as long as 
61 = (or #2 = 0), there are always zero modes trapped 
in the normal region of the junction involving iron pnic- 
tides with sign-changed s-wave pairing symmetry. How- 
ever, in a more realistic system with band structure such 
as our SNS junction it usually introduces finite effective 
mass for the band electrons, which destroys the validity 
of using linearized Eq. ([8]) (where the effective mass goes 
to infinity). As a result, the to-be-degenerate zero modes 
splits as we see in Fig. |3] (only E > shown). The same 
argument is also applicable when we further consider the 
effect brought by k y in our Hio(k y ). 
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